Transient time-domain resonances and the time scale for tunneling 
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Transient time-domain resonances found recently in time-dependent solutions to Schrodinger's 
equation are used to investigate the issue of the tunneling time in rectangular potential barriers. 
In general, a time frequency analysis shows that these transients have frequencies above the cutoff 
frequency associated with the barrier height, and hence correspond to non-tunneling processes. We 
find, however, a regime characterized by the barrier opacity, where the peak maximum tmax of the 
time-domain resonance corresponds to under-the-barrier tunneling. We argue that tmax represents 
the relevant tunneling time scale through the classically forbidden region. 
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Tunneling refers to the possibility that a particle tra- 
verses through a classically forbidden region. In the en- 
ergy domain the solution to Schrodinger's equation at 
a fixed energy E is a subject discussed in every quan- 
tum mechanics textbook. In the time domain, however, 
there are still aspects open to scrutiny. A problem that 
has remained controversial over the years is the tunnel- 
ing time problem, that may be stated by the question: 
How long does it take to a particle to traverse a classi- 
cally forbidden region? Different authors have proposed 
and defended different views in answering the above ques- 
tion Hiii. 

In recent work we have investigated the effect of the 
transient solutions to the time-dependent Schrodinger's 
equation for cutoff wave initial conditions (quantum shut- 
ter) on the tunnehng process [1,0,01- particular we 
found that just across the tunneling barrier, the proba- 
bility density as a function of time may exhibit a tran- 
sient steucture that we have named time-domain reso- 
nance. The peak value tmax of this structure represents 
the largest probability of finding the particle at the bar- 
rier width, L. More recently, in collaboration with Del- 
gado and Muga Q , we considered a time-frequency anal- 
ysis 0, to show the existence of under-the-barrier 
transients (forerunners)in very broad barriers. However, 
these occur along a finite region of the potential, and 
hence do not allow to characterize the time scale asso- 
ciated to the tunneling process through the full classi- 
cally forbidden region. On the other hand, in collabora- 
tion with Yamada jllj we have recently established the 
equivalence of our formulation with the notion of 'pas- 
sage time' in the real Feynman histories approach. The 
'passage time' yields the traversal time through a barrier 
region though it does not distinguish between processes 
from above or below the barrier height. The above con- 
siderations indicate that it is not clear under what con- 
ditions the time-domain resonances found at the barrier 
edge X — L correspond to genuine tunneling processes. 



The aim of this work is to show the existence of a 
regime, characterized by the opacity of the system, where 
the time-domain resonance maximum, t^axi corresponds 
to under-the-barrier tunneling. We argue that this time 
scale provides the tunneling time through the classically 
forbidden region. 

Our approach to the tunneling time problem is based 
on a model that deals with an explicit solution |5j to the 
time-dependent Schrodinger's equation 
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for an arbitrary potential V{x), defined in the region (0 < 
X < L), that vanishes outside that region. We consider 
the problem of the time evolution of a cutoff plane wave 
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following the instantaneous opening at t = of a quan- 
tum shutter at a; = 0. Along the tunneling region the 
solution reads, 
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where the 0±fc(x)'s refer to the stationary solutions of 
the problem and 4>n{x) — 2ikun{0)un{x) / (k^ ~^n)- Sim- 
ilarly, the solution '^'^{x,k;t) for the external or trans- 
mitted region {x > L), is given by [3, 

^%x,k;t) - TkMiyk) - T_fcM(y_fc) 

OO 

-^ J2 TnMivkJ, {x>L) (4) 

where the T±kS refer to the transmission amplitudes, 
and the factor r„ = 2ifcu„(0)M„(L) exp(— iA:„L)/(fc^ — 



2 



k^). In Eqs. lO and Q the coefficients Xn and T„ are 
given in terms of the resonant eigenfunctions {un{x)} 
with complex energy eigenvalues i?„ = h?k^/2m, with 
kn — o-n — ibn (a„,6n > 0). The resonant sums in Eqs. 
(O and ^ run over the full set of complex poles {fc„}. 



The M's are defined as [3, 
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where w is the complex error function jl2| defined as 
-iz), with arguments yq{x,t) = 



w{z) = exp(— z^)erfc(- 



e-^^/*(rn/2ht)^/'^[x - hqt/m], where q = ±fc, fc„. 

We shall begin exploring the issue of the tunneling time 
scale by first considering specific examples of rectangu- 
lar potential barriers (of height V and thickness L) to 
go then to results of a more general character. Let us 
recall the main features of a time-domain resonance |^. 
It follows from the initial condition, given by Eq. 
that initially there is no particle along the tunneling re- 
gion. Hence evaluating the probability density at the 
barrier width x = L as time evolves from zero, yields a 
distribution of characteristic times associated to the tun- 
neling process. In Fig. ^we plot the probability density 
^''^P (solid line), normalized to the transmission coef- 
ficient |Tfcp, as a function of time t, correspondin g t o 
a potential barrier system with typical parameters [13|: 
V = 0.3 eV, incidence energy E = h'^P/2m = 0.001 eV, 
and effective mass for the electron m = O.OGTm^, for a 
fixed value of position, x ^ L = 4.0 nm. We can clearly 
appreciate a structure peaked at the value tmax — 5.17 
fs. This is the so called time-domain resonance. 

In Fig. 121 we plot tmax (full dot) for the same parame- 
ters as in the previous figure except for the barrier width 
L, that we vary. Here we can clearly observe the exis- 
tence of a basin along a range of values of the barrier 
width. We can also appreciate that if L is further in- 
creased, tmax starts to grow linearly with L. Such a linear 
regime occurs at large barrier widths. We have discussed 
elsewhere, that this last situation refers to non-tunneling 
processes We have also pointed out in Ref. 6] that 
for small values of the barrier width L, the basin exhib- 
ited by trnax, is the result of a subtle interplay between 
tunneling and (non-tunneling) top-barrier resonant pro- 
cesses. 

In what follows we shall investigate under what condi- 
tions the time scales associated to the basin, are in fact 
related to a genuine tunneling process. 

We use the fact that the initial cutoff wave possesses 
a distribution of momentum components in k-space, and 
hence also of frequency components. As time evolves and 
the wave interacts with the potential, these frequency 
components manifest themselves in the time evolution 
of the probability density. The frequency content of the 
time-domain resonance can be investigated by perform- 
ing a time-frequency analysis. We do this by computing 
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FIG. 1: Time evolution of l*""]^ (solid line), at a; = L = 4.0 
nm. A full square indicates the position of the maximum of 
the time-domain resonance, at tmax = 5.17 fs. 
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FIG. 2: Maximum of the time-domain resonance tmax (full 
dot) as a function of the barrier width L, for an incidence 
energy E — .001 eV. In this case the barrier height isV = 0.3 
eV. See text. 



the local average frequency uJav ^3 > 

c^au = -Im 
and the instantaneous bandwidth a jinl |. 
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where s = i,e refers, respectively, to the internal and ex- 
ternal solutions. To exemplify this, we choose the case 
L — 4.0 nm depicted in Fig. El which is located around 
the minimum of the basin. In Fig. (31 we plot the rela- 
tive average local frequency (relative frequency for short) 
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FIG. 3: Relative average local frequency ujav/i^v (solid 
line) for the case depicted in Fig. Q The cutoff-frequency 
uiav/i^v = 1 (dashed line) is included for comparison. In the 
inset we plot the instantaneous bandwidth a of the spectro- 
gram depicted in the main graph. Notice that the frequency 
deviations at tmax are exactly zero, i.e., a{tmax) = 0. In all 
cases a full square indicates the position of tmax- 



^av/^v, where ujv = V/h is the cutoff frequency, along 
the relevant time interval, discussed in Fig. We can 
appreciate that in the vicinity of the maximum of the 
time-domain resonance, t„iax, the probabiUty density is 
composed entirely by under-the-barrier frequency com- 
ponents i.e ojav/^^v < 1- This also occurs at the ex- 
act value tmax, also indicated in the figure by a solid 
square. In the inset of Fig. |2|we plot the instantaneous 
bandwidth a of the spectrogram. Notice the absence of 
a frequency dispersion around the maximum tmax, *-e., 
<^{tmax) = 0. The above result indicates that in this case 
the peak of the time-domain resonance, and the values 
close to it, refer to a tunneling event. We have found, 
however, that this is not a general situation. For in- 
stance, for values of L outside the basin, i.e., along the 
linear regime in Fig. [2Jl, the average frequency related 
to the corresponding tmax is above the cutoff frequency 
and hence refers to non-tunneling processes. As we shall 
present below this is more appropriately discussed by us- 
ing the notion of the opacity of the system. In Fig. 21 we 
plot the relative frequency ujav/^v associated to differ- 
ent values of the maximum tmax , measured at different 
positions along both the internal and external regions of 
a potential barrier with parameters: V — 0.3 eV, and 
L = 4.13 nm. In this case we choose the following values 
of the incidence energy: E = 0.001 eV (solid dot), and 
E = 0.01 eV (hollow dot). In the inset of Fig. ^we show, 
for the particular case oi E = 0.01 eV, the values of tmax 
(solid square) at the different values of position consid- 
ered in the main graph. As can be clearly appreciated in 
that figure, the tunneling process along the whole inter- 
nal region is governed by under-barrier-frequency com- 



FIG. 4: Relative frequency u>av/i^v of the maximum of the 
time-domain resonance, as a function of position, measured in 
units of the barrier width L. The parameters are given in the 
text. Two incidence energies are considered: E — 0.001 eV 
(solid dot), and E = 0.01 eV (hollow dot). In both cases, the 
relative frequency ujav/^^v along the internal region is below 
the cutoff- frequency LOavj^v ~ 1 (dashed line). The behavior 
of tmax as a function of position is illustrated in the inset for 
the case with E = 0.01 eV. The position of the barrier edge, 
a; = L, is indicated by a dotted line in both figures. 



ponents, i.e., ujavl^v < 1- We can see in Fig. ^ that we 
can still observe frequency components below the cutoff- 
frequency, LOy for distances up to a: ~ 2L along the ex- 
ternal region. As xjh increases further, LOav/^v > 1- 
This behavior indicates the prevalence of non-tunneling 
components in the behavior of the probability density . 

We have found that the regime which corresponds to 
under-the-barrier time-domain resonances at the barrier 
edge X ^ L, may be described more generally by referring 
to the opacity a of the system, defined as. 



[2mV] 



1/2 



-L, 
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and by the dimensionless parameter u, the ratio between 
the potential barrier height and the incidence energy. 



V 
E' 



(9) 



To characterize this tunneling regime, we use the fact 
that all systems sharing the same parameters a and u 
yield the same relative frequency, cuav/^v- This regular- 
ity arises from a simple rescaling property of the time- 
dependent Schrodinger's equation and the corresponding 
initial condition. By feeding the dimensionless variables 
X = x/L, and T = uiyt in Eqs. Q and (O, we obtain. 
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FIG. 5: Relative frequency Uav/i^v measured at the barrier 
edge a; = L, as a function of the opacity a. Here we considered 
a barrier height V — 0.3 eV, and the parameters: u = 300 
(solid dot), II = 10 (solid triangle), and u — 5 (solid square). 
Note that for values of the opacity smaller than a ~ 3.3, the 
relative frequencies for all values of u are below the cutoff- 
frequency Uav/i^v = 1 (dashed line). See text. 



with the initial condition 



X{X,T^0) = 



-iaX/ -Ju 



0, 



X < 0, 
X > 0, 



(11) 



where x(^, is the rescaled time dependent-solution. 
From Eqs. Eqs. I|10|l and it is clear that the time- 
dependent solution must depend only on the parameters 
a and u. That is, for a fixed value of a, all the systems 
with the same parameter it, yield the same x{X^T^. As 
a consequence of the above considerations we can write 
the relative frequency as, 
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and use it to characterize the regime associated with 
under-the-barrier frequency components. In particular, 
we are interested in defining the range of values of a 
where the relative frequencies associated to the time- 
domain resonance are below the cutoff- frequency . In 
Ref. it is found that for opacities less than a criti- 
cal value no time domain resonances occur. We denote 
it by amin and it has the value a„im = 2.065. In Fig. 
iniwe plot the relative frequency tOav/^^v as function of 
the opacity a, for three different values of the param- 
eter u, u = 5 (solid dot), u = 10 (solid triangle), and 
u — 300 (solid square). In this case we have chosen a 
value of V = 0.3 eV in the calculation. Although in Fig. 
0we have considered values of the parameter u such that 
5 < u < 300, the cases corresponding to very large values 
of M (m — > oo) (not shown here), almost overlap with the 
case u = 300. Thus, for very large values of u we find 



a maximum value for the opacity amax ~ 3.3. Conse- 
quently one may define an opacity "window" , in the range 
of values 2.065 < a < 3.3, where the relative frequencies 
are always below the cutoff- frequency ujav/uJv — 1, irre- 
spective of the value of the parameter m, namely of the 
value of the incidence energy. Note that the above nu- 
merical values refer to the effective mass m — 0.067me 
and clearly will be modified for other values of the effec- 
tive mass. 

It is of relevance to point out that along the basin re- 
gion the time scale given by the peak maximum tmax 
differs in a essential way from both the semi-classical 
Biittiker-Landauer and Biittiker traversal times, which 
in addition to exhibit always a linear dependence with 
L, refer to over-the-barrier processes Also tmax 

represents a completely different notion than the phase- 
time, which corresponds to a long-time asymptotic no- 
tion representing a global effect of the potential on the 
Schrodinger's solution, as discussed in Refs. 0,0- 

To conclude we remark that the analytical solution to 
the time-dependent Schrodinger's equation with quan- 
tum shutter initial conditions applies in general to arbi- 
trary potentials, provided they vanish beyond a distance, 
and can also be extended to deal with finite cutoff pulses 
as discussed in Ref. 7]. The quantum shutter setup 
provides a consistent procedure to obtain the timneling 
time: initially there is no particle along the tunneling re- 
gion and as time evolves the transient peaked structure 
exhibited by the probability density at the barrier width 
provides the relevant time scale for tunneling. This oc- 
curs within a range of values of the opacity of the system 
and is independent of the incidence energy. It is worth 
noticing that the values of a within the opacity "window" 
may be obtained using typical parameters of semiconduc- 
tor heterostructures [l3|. Also one should stress that at 
the peak maximum, the time-domain resonance is gov- 
erned by a single frequency. That is, the system acts as a 
frequency filter. To test our results experimentally would 
require to consider the detection of tunneling particles in 
time domain at distances close to the interaction region. 
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